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Introducción 

La  imagen  adquirida  por  un  sensor  a  bordo  de  un 
satélite  presenta  diferencias  con  respecto  a  la  imagen 
real  tanto  en  la  situación  de  sus  puntos  como  en  el 
nivel  de  brillo  de  cada  uno.  En  consecuencia,  la  geo¬ 
metría  de  la  imagen  no  refleja  con  fidelidad  absoluta 
la  de  la  superficie,  ni  tampoco  pueden  deducirse  sus 
propiedades  físicas  a  partir  de  la  lectura  directa  de 
los  niveles  digitales.  Es  necesario  aplicar  unas  deter¬ 
minadas  operaciones  de  corrección,  que  están  enca¬ 
minadas  a  minimizar  estas  alteraciones  en  la  medida 
de  lo  posible.  Algunas  de  estas  tareas  son  aplicadas 
por  las  mismas  estaciones  de  recepción,  las  cuales 
ofrecen  distintos  niveles  de  corrección.  El  estudio  de¬ 
tallado  de  la  imagen,  los  estudios  multitemporales  o 
la  fusión  con  otras  fuentes  cartográficas  requieren 


correcciones  minuciosas,  que  normalmente  son  reali¬ 
zadas  por  el  analista. 

La  georreferenciación  es  el  proceso  mediante  el  cual 
se  dota  a  la  imagen  de  coordenadas  cartográficas,  en 
el  sistema  de  proyección  que  se  elija,  tarea  que  suele 
abordarse  simultáneamente  con  la  corrección  geomé¬ 
trica. 

Este  capítulo  aborda  las  operaciones  más  habituales 
de  corrección  geométrica  en  el  pretratamiento  de  las 
imágenes  de  satélite,  operación  necesaria  antes  de 
proceder  al  estudio  en  sí  de  la  información  que  pro¬ 
porcionan. 


Métodos  polínómicos  de  corrección 


La  corrección  geométrica  de  las  imágenes  requiere 
deducir  una  función  de  transformación  para  cada  va¬ 
riable  planimétrica  que  convierta  la  posición  inicial 
{x,y)  de  cada  celda  en  coordenadas  corregidas  (n,v): 


v  =  gl(x,y)  J 

Que  también  se  pueden  poner  en  función  de  las  coor¬ 
denadas  imagen  en  columna  y  línea  de  cada  celda: 


u  =  f2(c,l) 

v  =  g2{c,l) 

Dado  que  las  distorsiones  geométricas  sufridas  por  la 
imagen  obedecen  a  numeras  causas,  algunas  deriva¬ 
das  del  efecto  proyectivo,  de  la  esfericidad  terrestre, 
de  la  rotación  del  planeta,  etc.  cuyo  análisis  puede 
abordarse  analíticamente,  otras  tantas  son  aleatorias 
y  vienen  a  engrosar  y  complicar  el  efecto  resultante. 
En  definitiva,  la  distorsión  geométrica  estaría  defi¬ 
nida  por  una  función  complicada  y  desconocida. 
Cuando  no  se  conoce  explícitamente  una  función  di- 
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ficultosa  de  estas  características,  se  le  suele  caracte¬ 
rizar  matemáticamente  mediante  su  desarrollo  en  se¬ 
rie,  de  cuyo  polinomio  se  toman  sólo  los  términos  ne¬ 
cesarios,  sabiendo  que  un  cierto  número  de  ellos  pro¬ 
porcionan  una  explicación  razonablemente  suficiente 
del  comportamiento  de  la  función.  En  consecuencia, 
a  efectos  prácticos,  se  puede  sustituir  la  función  ori¬ 
ginal  por  dicho  polinomio. 

Los  métodos  más  sencillos  de  corrección  geométrica 
están  basados  en  la  adopción  de  funciones  polinómi- 
cas  que  transforman  las  coordenadas  de  un  conjunto 
de  puntos  de  control  en  coordenadas  corregidas.  Una 
vez  encontrados  los  polinomios,  se  aplican  a  todas  las 
celdas  de  la  imagen  para  forzar  su  deformación  hasta 
conseguir  una  imagen  corregida  geométricamente. 

El  procedimiento  sigue  las  etapas  siguientes: 

—i  Selección  de  los  puntos  de  control. 

-i  Cálculo  de  los  coeficientes  de  los  polinomios 
de  ajuste. 

—i  Transferencia  de  las  celdas  a  sus  nuevas  posi¬ 
ciones. 

—i  Remuestreo. 

Selección  de  los  puntos  de  control 

Los  puntos  de  control  son  puntos  de  los  que  se  co¬ 
nocen  sus  coordenadas  reales  (coordenadas-terreno) 
y  las  coordenadas  matriciales  de  sus  homólogos  en  la 
imagen  (coordenadas-imagen). 

La  correcta  selección  de  los  puntos  de  control  y  la 
determinación  de  sus  coordenadas  reales  e  imagen  es 
determinante  en  el  resultado  final  del  proceso.  Los 
puntos  de  control  deben  ser  elegidos  en  accidentes 
del  terreno  permanentes,  como  cruces  de  caminos, 
por  ejemplo,  u  otros  hitos  visibles.  No  son  válidos  las 
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formas  definidas  por  sombras  en  la  imagen,  como 
tampoco  lo  son  las  derivadas  de  accidentes  hidrológi¬ 
cos,  puesto  que  su  posición  no  es  seguro  que  sea  es¬ 
table. 

La  calidad  del  ajuste  depende  de  tres  factores: 

1.  Número  de  puntos.  El  mínimo  número  nece¬ 
sario  para  resolver  las  ecuaciones  es  de  3  pun¬ 
tos,  si  los  polinomios  son  de  primer  grado,  6, 
si  son  de  segundo  grado  y  10  si  son  de  tercer 
grado.  Aunque  se  aconseja  tomar  no  menos 
del  doble  de  número  matemáticamente  nece¬ 
sario  para  el  ajuste,  siempre  es  preferible  to¬ 
mar  un  número  de  puntos  considerablemente 
mayor. 

2.  Localización  de  los  puntos.  Es  imprescindible 
un  gran  rigor  en  la  localización  geométrica 
exacta,  tanto  en  la  imagen  como  en  el  terreno 
o  en  la  fuente  cartográfica  que  lo  sustituya. 

3.  Distribución  de  los  puntos  en  la  imagen.  Es 
deseable  una  distribución  uniforme.  Por  otra 
parte,  siempre  es  preferible  la  interpolación  a 
la  extrapolación.  Con  ello  se  quiere  decir  que 
si  se  toman  puntos  periféricos  en  la  imagen 
puede  asegurarse  con  algún  margen  de  segu¬ 
ridad  que  el  ajuste  interior  pueda  ser  tolera¬ 
ble.  Sin  embargo,  la  inversa  no  es  cierta:  por 
muchos  puntos  que  se  tomen  en  una  zona 
central  de  la  imagen,  nunca  se  tiene  la  segu¬ 
ridad  de  que  las  zonas  periféricas  de  la  ima¬ 
gen  ajusten  correctamente,  especialmente  si 
se  utilizan  polinomios  de  segundo  o  tercer 
grado. 

Las  funciones  que  transforman  las  coordenadas  ac¬ 
tuales  ( x,y )  de  los  puntos  de  control  en  coordenadas 


corregidas  ( u,v )  del  mapa  de  referencia  son  polino¬ 
mios,  cuya  forma  genérica  es: 


n  n-p 

P=°  9=°  /0x 

f  (2) 

n  n-p 

P= 0  9=0 

donde  n  representa  el  grado  del  polinomio.  Los  poli¬ 
nomios  de  primer  grado  permiten  rotar,  trasladar  la 
imagen  y  cambiar  su  escala.  Los  de  segundo  incluyen 
ciertas  deformaciones  residuales  aleatorias.  Los  de 
tercer  grado  comportan  una  corrección  más  depu¬ 
rada.  No  es  habitual  utilizar  polinomios  de  grado  su¬ 
perior,  aunque  en  determinadas  aplicaciones  de  co¬ 
rrección  asistida  por  ordenador  se  utilizan  polino¬ 
mios  de  cuasi-cuarto  grado  y  cuasi-sexto  grado: 


u=Y¿HapqxPyq 

p= 0  q= 0 

v=Tj'éjbpqxPy9 

p= 0  q= 0 


(3) 


En  cualquiera  de  los  casos,  Los  coeficientes  apq  y  bpq 
se  determinan  por  ajuste  mínimo-cuadrático  a  partir 
de  las  coordenadas  primitivas  y  corregidas  de  los 
puntos  de  control.  El  número  mínimo  de  puntos  de 
control  que  deben  seleccionarse  para  encontrar  una 
solución  depende  del  grado  del  polinomio  y  sigue  la 
regla: 


m  = 


(n  +  l)(n  +  2) 


(4) 


Donde  m  es  el  mínimo  número  de  puntos  de  control 
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y  n  el  grado  del  polinomio  adoptado  para  el  ajuste. 

Siguiendo  el  formato  (2),  El  caso  más  sencillo  es  el 
de  ajuste  lineal,  en  el  que  las  coordenadas  corregidas 
(w,v)  se  encuentran  en  función  de  las  coordenadas 
{x,y)  de  los  puntos  de  control  de  la  imagen  sin  corre¬ 
gir: 


u  =  a0  +alx  +  a2y 
v  =  b0  +blx  +  b2y 


(5) 


Al  ser  3  las  incógnitas  para  cada  una  de  las  dos  varia¬ 
bles  ao,  <21,  ai  y  bo,  b\,  bi,  su  resolución  requiere  al 
menos  de  un  sistema  de  3  ecuaciones.  Se  necesitarían 
por  lo  tanto  un  mínimo  de  3  puntos  de  control  para 
este  ajuste  y  ello  podría  proporcionar  una  solución 
matemática,  aunque  no  necesariamente  correcta 
desde  el  punto  de  vista  cartográfico.  Siempre  que  el 
sistema  de  ecuaciones  sea  compatible  se  puede  dar 
con  una  solución,  incluso  habiendo  tomado  alguna 
coordenada  erróneamente.  Sin  embargo,  el  ajuste  mí¬ 
nimo  cuadrático  se  basa  en  la  sobreabundancia  de 
datos  para  encontrar  una  solución  que  satisfaga  ra¬ 
zonablemente  al  conjunto  de  todos  ellos.  Por  esta  ra¬ 
zón  se  procederá  siempre  a  seleccionar  un  número 
no  menor  del  doble  del  mínimo  matemáticamente  ne¬ 
cesario,  aunque  es  deseable  un  número  considerable¬ 
mente  superior. 

Para  el  ajuste  mediante  polinomios  de  segundo 
grado,  las  ecuaciones  a  plantear  serán: 

u  =  a0  +  a{x  +  a2y  +  arxy  +  a4x2  +a5y 2 
v  =  b0  +  bxx  +  b2y  +  b2xy  +  b4x2  +b5y 2 

Para  la  resolución  de  los  coeficientes  se  necesitará 
un  sistema  de  6  ecuaciones  con  6  incógnitas  para 
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cada  una  de  las  variables.  Por  ello,  el  mínimo  número 
de  puntos  de  control  a  utilizar  en  un  ajuste  de  se¬ 
gundo  grado  deberá  ser  de  6,  pero  en  la  práctica  no 
se  tomarán  menos  de  12  puntos. 

Un  ajuste  de  tercer  grado  exigiría  plantear  el  si¬ 
guiente  sistema: 

u  =  a0  +  axx  +  a2y  +  a2xy  +  ci4x2  +  a5y2  +  a6xy2  + 
+a1x2y+  a&x3  +  a9y3 

7  f  (7) 

v  =  b0  +  bxx  +  b2y  +  b2xy  +  b4x2  +  b5y2  +  b6xy~  + 

+b1x2  y  +  \x 3  +  bg  y3 

para  cuya  resolución  se  requiere  un  mínimo  de  10 
puntos,  aunque  en  la  práctica  se  tomarán  al  menos 
20  puntos.  En  correcciones  de  imágenes  completas 
se  utilizarán  por  encima  de  50  o  60  puntos,  bien  dis¬ 
tribuidos  por  la  escena. 

En  el  caso  de  utilizar  un  polinomio  de  cuasi-cuarto 
grado,  el  sistema  de  ecuaciones  sería: 


u  =  a0  +  axx  +  a2y  +  a3xy  +  aAx2  +  a5y2  +  a6xy2  + 

+a7x2  y  +  axx2  y2 

,  (8) 

v  =  b0  +  bxx  +  b2y  +  b3xy  +  b4x  +  b5y2  +  b6xy2  + 

+b7x2y  +  b&x2y2 

Como  se  ve,  con  este  formato  se  requieren  al  menos 
9  puntos  de  control,  aunque  se  seguirán  tomando 
como  poco  el  doble. 


Ajuste  mínimo-cuadrático 

Como  ha  quedado  dicho,  para  resolver  los  coeficien¬ 
tes  de  un  polinomio  de  grado  n  no  sirve  plantear  un 
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sistema  de  tantas  ecuaciones  como  coeficientes.  Si  el 
sistema  es  compatible,  siempre  arrojaría  una  solu¬ 
ción,  con  independencia  de  la  correcta  o  incorrecta 
selección  de  los  puntos  o  de  la  introducción  de  sus 
coordenadas  reales.  Para  asegurar  el  cálculo  razona¬ 
blemente  correcto  de  los  coeficientes,  se  recurre  a 
ajuste  mínimo-cuadrático. 

En  el  ajuste  por  mínimos  cuadrados  se  seleccionan 
un  número  m  de  puntos  considerablemente  mayor 
del  necesario  y  se  plantean  con  ellos  un  sistema  de 
m  ecuaciones  con  n  incógnitas. 

Para  m  puntos,  el  sistema  de  ecuaciones  mediante 
polinomios  de  segundo  grado,  sería,  siguiendo  (6): 
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o  sea: 


(9) 


A-B„=C  1 
A-BV=CVJ 


(10) 
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A  B  =  C 


(’  c  tiene  igual  numero  de  filas  que  a 
e  Igual  número  de  columnas  que  B 


Xc 

A  rw  es  cwij'sdo  y  por  lo 


Procedimiento  para  obtener  el  vector  B. 


donde  A  es  la  matriz  de  coordenadas  imagen  de  los 
puntos,  B  el  vector  de  coeficientes  y  C  el  vector  de 
coordenadas  terreno. 

Como  A  no  es  cuadrada,  no  es  invertible,  por  lo  que 
hay  que  recurrir  a  un  pequeño  artificio  matemático 
para  poder  despejar  B,  que  es  lo  que  interesa.  Multi¬ 
plicando  ambos  miembros  por  la  matriz  traspuesta  de 
A  se  consigue  en  el  primer  miembro  de  cada  una  de 
las  expresiones  (10)  una  matriz  que  sí  puede  ser  in¬ 
vertida  y  pasada  al  segundo  miembro,  quedando: 


b„=(a'a)‘V  c, 
B,=(A'A)JA'  c, 

que  es  la  solución  buscada. 


(11) 


El  enfoque  del  problema  desde  la  perspectiva  de  los 
mínimos  cuadrados  conduce  al  mismo  resultado.  Se 
conocen  las  variables  x  e  y,  que  son  las  coordenadas 
imagen  de  los  puntos  de  control.  Las  incógnitas,  por 
el  contrario,  son  los  coeficientes  de  los  polinomios  a¡ 
y  bu  Tampoco  se  conocen  las  coordenadas  u  y  p  de 
cada  punto  corregido,  puesto  que  todavía  se  desco¬ 
nocen  los  coeficientes  que  conforman  cada  polino¬ 
mio. 


Xx  =  a0  +  axxx  +  a2yx 
X2  =  a0  +axx  2  +a2y2 
X3  =  a0  +axx 3  +  a2y3 
X4  =  a0  +  alx4+a2y4 

que  en  forma  matricial  sería: 


O 

1 

1 

1 


T2 

y4J 


ir,  \ 


a , 


\a2  J 


*2 

V*47 


(12) 


(13) 


Pero  en  realidad  sucede  que: 

a0  +  a]x[  +a2yx  =  a, 
a0  +  a]x2  +a2y2  =  u2 
a0  +axx3  +a2y3  =  u3 
a0  +axx4  +  a2y4  =  u4 

que  difiere  de  la  posición  teórica  X¡. 
Calculando  el  error  cometido  para  cada  punto: 


(14) 


Para  formalizar  el  método,  de  cada  punto  se  determi¬ 
nan  sus  coordenadas  reales  X  e  Y,  que  son  las  coor¬ 
denadas  terreno.  La  diferencia  entre  la  posición  (u,v) 
y  (X,Y)  es  el  error  del  ajuste.  Se  pretende  que  la  suma 
de  los  cuadrados  de  todos  ellos  sea  mínima. 

Supóngase  un  ajuste  de  primer  grado,  que  sólo  exi¬ 
giría  3  puntos  y  4  puntos  de  control.  Asumiendo  que 
el  ajuste  fuera  perfecto: 


s,  =  tq  -  Xí  =  a0  +  axx j  +  a2yx  -  Xx 
e-,  =u1-X1  =  a0  +axx2  +a2y2  -X0 

> 

s 3  =  w3  -  =  a0  +  axx3  +a2y3-  X3 

s4  =u4-X4  =a0  +  axx4  +a2y4-X4 

y  de  igual  forma  para  la  Y. 

Se  suman  los  cuadrados  de  todos  los  errores: 
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°  =  Z 8<2  =  Z K  +  «i*;  +  «2  Y  -  *,■  )2  ( 16) 

1=1 

y  se  impone  la  condición  de  mínimo  igualando  a  cero 
la  primera  derivada: 


50 

da0 

50 

dax 

50 

da2 


4 

2X(ao+ai  Xi+a2yi~xi)  =° 

i=l 

4 

2Z  K + ai xi  +  a2 y ¡  -  xi )  x,  =  o  > 

i= l 
4 

2X  K  +  aiXi  +  =  O 

í=l 


Operando: 


4a0  +  alYJxi+a2Yjyi=YjXi 
a0  ]T  x,.  +a,X  xf  +  «2  Z  Y  =  Z  > 
ao  Z  ><■ + «i  Z  a  y  i  +  a2  Z  y?  =  Z  y . 


(17) 


(18) 


Y  poniendo  la  expresión  en  forma  matricial: 


f  4 

2>, 

2>.  j 

(a  > 
Mo 

rz*,  i 

2>< 

2><?, 

ai 

= 

z** 

2>,y, 

2>,2 , 

\a2  y 

A‘A 


que  en  realidad  es: 

A'AB  =  A‘C  (20) 

es  decir: 

B  =  (a‘a)'1A‘C  (21) 

expresión  análoga  a  (11). 


Otros  métodos  de  corrección 


Aunque  los  más  utilizados  son  los  métodos  de  correc¬ 
ción  polinómica,  existen  otras  formas  de  asemejar  la 
imagen  a  un  mapa,  y  el  analista  deberá  valorar  si  con¬ 
vienen  emplearlos.  A  continuación  se  exponen  breve¬ 
mente  algunos  de  ellos. 

Transformación  afín 


coordenadas  de  dos  puntos  en  el  sistema  de  referen¬ 
cia  de  la  imagen  y  en  el  del  mapa,  es  decir,  cuatro 
coordenadas  en  total.  Se  conoce  también  como  trans¬ 
formación  de  «4  parámetros»,  que  son:  Eo,  No,/,  £,  es 
decir,  las  dos  traslaciones  este  y  norte  del  origen  de 
coordenadas,  el  factor  de  escala  y  el  ángulo  de  rota¬ 
ción,  respectivamente. 


Consiste  en  la  combinación  de  una  rotación,  una  tras¬ 
lación  y  un  cambio  de  escala.  Requiere  conocer  las 


vvy 


í  p  A 


+  / 


^coss 

vsens 


-sen 
eos  s  j 


(22) 
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Transformación  afín. 


donde  u  y  v  son  las  coordenadas  corregidas  de  un 
determinado  punto  y  x  e  y  son  sus  coordenadas  ori¬ 
ginales. 

Los  cuatro  parámetros  se  calculan  mediante  un 
ajuste  por  mínimos  cuadrados  de  un  conjunto  de 
ecuaciones  bastante  superior  a  cuatro. 

Esta  transformación  es  demasiado  simple  para  corre¬ 
gir  las  deformaciones  habituales  de  una  imagen  de 
satélite.  La  razón  se  debe  al  hecho  de  que  el  algo¬ 
ritmo  simplemente  escala  el  mapa  por  un  factor,  pero 
no  modifica  la  matriz  de  celdas.  Además,  el  factor  de 
escala  es  único  para  toda  la  escena,  lo  que  exigiría 
que  la  imagen  hubiera  sido  deformada  por  igual  en 
toda  su  extensión  y  en  todas  las  direcciones. 


Control  de  calidad  del  ajuste 

Las  funciones  polinómicas  de  corrección  definen  lo 
mejor  posible  los  cambios  de  posición  que  han  de  su¬ 
frir  los  puntos  de  control.  Estos  cambios  no  obedecen 
en  términos  generales  a  transformaciones  matemáti¬ 
cas  simples,  debido  a  la  aleatoriedad  de  las  distorsio¬ 
nes  sufridas  por  la  imagen. 

Las  funciones  encontradas  son  aquéllas  que  mejor 
satisfacen  las  posiciones  de  todos  los  puntos  de  con¬ 
trol,  estableciendo  un  compromiso  entre  ellos,  de  tal 
modo  que  la  suma  de  los  cuadrados  de  las  desviacio¬ 
nes  sea  mínima.  El  indicador  más  utilizado  para  me¬ 
dir  la  calidad  del  ajuste  es  el  error  medio  cuadrático 
RMS  (Root  Mean  Square ): 
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Rubber-sheeting 

La  transformación  denominada  rubber-sheeting 
toma  su  nombre  de  su  similitud  con  la  adaptación  de 
una  lámina  elástica  a  un  mapa  base  indeformable.  La 
imagen  de  satélite  estaría  impresa  sobre  la  superficie 
de  un  globo  hinchable  de  caucho,  con  multitud  de 
deformaciones,  y  se  trataría  de  hacer  coincidir  un 
conjunto  de  puntos  esenciales  de  la  imagen  -de  con¬ 
trol—,  sujetándolos  con  alfileres  o  chinchetas  al  pro¬ 
pio  mapa.  Realmente  se  trata  de  un  modelo  polino¬ 
mio)  de  corrección  basado  en,  al  menos,  cuatro  pun¬ 
tos: 

u  =  a0  +  axx  +  a2y  +  a3xy 
v  =  b0  +  bvx  +  b2y  +  b3xy 


RMS  = 


ZR«-x,)2+(v,-r)2 


1=1 


(24) 


en  cuya  expresión,  (X,  Yi)  son  las  coordenadas  reales 
de  los  puntos  de  control,  (ui,  vi)  son  las  coordenadas 
estimadas  para  los  puntos  de  control  a  partir  de  la 
aplicación  de  las  funciones  de  ajuste  y  n  es  el  número 
de  puntos  de  control  utilizados. 

En  la  práctica  se  admite  un  RMS  máximo  de  2/3  del 
tamaño  de  la  celda  para  cada  coordenada  X  e  Y,  pero 
se  pueden  admitir  desviaciones  de  cada  punto  equi¬ 
valentes  al  tamaño  de  la  celda  tanto  en  X  como  en  Y, 
ya  que  es  el  límite  de  la  resolución  espacial  de  la  ima¬ 
gen. 
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Ortorrectificación 


El  relieve  del  terreno  tiene  un  carácter  caprichoso  y 
su  mayor  o  menor  aleatoriedad  no  responde  a  nin¬ 
guna  función  sistemática.  Como  consecuencia,  las 
distorsiones  geométricas  que  introduce  en  la  imagen 
persisten  tras  la  corrección  geométrica  polinómica. 
Para  eliminar  el  desplazamiento  debido  al  relieve  es 
necesario  someter  la  imagen  a  un  proceso  de  rectifi¬ 
cación  diferencial  punto  a  punto,  si  bien,  cuando  el 
desplazamiento  de  cualquier  celda  sea  inferior  a  la 
resolución  espacial  de  la  imagen,  podría  prescindirse 
de  esta  corrección. 

En  una  proyección  cónica  el  desplazamiento  debido 
al  relieve  que  sufre  un  elemento  de  una  altura  Ah 
sobre  el  plano  de  comparación  situado  a  una  cierta 
distancia  del  nadir  de  la  imagen  sigue  la  relación: 


A  r  Ah 
r  h 


(25) 


donde  h  es  la  altura  del  instrumento  de  adquisición, 
Ah,  como  se  ha  dicho,  es  la  altura  del  objeto,  r  es  la 
distancia  del  extremo  proyectado  medida  desde  el 
punto  principal  de  la  imagen  y  A r  es  el  desplaza¬ 
miento  debido  al  relieve. 


La  corrección  implica  conocer  el  desplazamiento  de¬ 
bido  al  relieve  despejando  de  (25): 


*  A/i 
A  r  =  r  — 
h 


(26) 


La  ortorrectificación  es  el  proceso  mediante  el  cual 


se  corrige  el  desplazamiento  de  cada  una  de  las  cel¬ 
das  de  la  imagen,  en  función  de  su  altura  relativa  y 
su  distancia  a  la  traza  del  satélite. 

La  expresión  (26)  puede  particularizarse  para  cada 
una  de  las  coordenadas  de  la  siguiente  forma: 


'  A h'] 

X  =  x  -  Ax  =  x 

v  h 

f  Ah 

Y  =  y- Ay  =  y 

1 - 

V  h  v 

(27) 


Pero  las  imágenes  de  satélite,  que  no  son  adquiridas  Geometría  de  la  proyección  cónica, 

por  fotogramas  completos  sino  por  líneas  o,  a  lo 
sumo,  por  pequeños  grupos  de  ellas,  no  se  adaptan  a 
una  proyección  cónica  sino  más  bien  a  una  de  tipo 
cilindrico  donde  las  distorsiones  causadas  por  la  to¬ 
pografía  adquieren  direcciones  perpendiculares  a  la 
traza  del  satélite  y  son  proporcionales  a  la  distancia 
del  punto  hasta  dicha  traza.  En  consecuencia,  puede 
prescindirse  de  una  de  las  dos  componentes.  Si  se 
considera  que  la  traza  mantiene  la  dirección  Y,  el  des¬ 
plazamiento  a  corregir  solamente  será: 


f 

X  =  x-Ax  =  x  1 

V 


A h' 
~h  ; 


(28) 


En  la  práctica,  el  valor  Ah  se  toma  de  un  modelo  di¬ 
gital  del  terreno  (MDT). 
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Imagen  original 


Remuestreo 


Imágenes  transformadas 


a  Igual  escala  con  reducción  con  aumento 


Necesidad  del  remuestreo. 


imagen  transformada 


imagen  corregida 


7 

i  / 

\ 

Tras  ser  sometidas  al  cambio  de  coordenadas  me¬ 
diante  la  transformación  polinómica  elegida  y,  en  su 
caso,  corregido  también  el  efecto  del  relieve,  las  posi¬ 
ciones  de  las  celdas  ocuparán  valores  reales,  que  en 
general  no  serán  números  enteros.  Pero  para  que  la 
imagen  pueda  ser  procesada,  se  requiere  que  cada 
celda  corregida  ocupe  una  posición  entera  de  co¬ 
lumna  y  línea. 

En  el  tratamiento  digital  de  imágenes  se  denomina 
remuestreo  ( resampling )  al  conjunto  de  operaciones 
necesarias  para  asignar  nuevos  valores  a  las  celdas 
cuando  se  produzca  alguna  circunstancia  como  con¬ 
secuencia  de  la  cual  sus  posiciones  no  sean  números 
enteros.  Un  ejemplo  de  tales  circunstancias  pueden 
ser  el  cambio  de  resolución  espacial  de  una  imagen, 
mediante  el  que  se  cambia  el  tamaño  final  de  la  celda, 
que  no  tiene  que  ser  necesariamente  múltiplo  ni  sub¬ 
múltiplo  del  tamaño  original.  En  lo  atañe  a  la  correc¬ 
ción  geométrica,  el  resultado  de  la  transformación  po¬ 
linómica  obliga  a  realizar  un  remuestreo  final  antes 
de  hacer  utilizable  la  imagen.  Mediante  el  remues¬ 
treo,  la  imagen  transformada  se  adapta  a  la  rejilla  del 
sistema. 

Al  tener  que  situar  el  centro  de  la  celda  en  una  posi¬ 
ción  intermedia  a  la  que  proporciona  la  transforma¬ 
ción,  los  nuevos  ND,  en  coordenadas  de  la  imagen  de 
salida,  se  recalculan  en  función  de  los  de  su  entorno, 
existiendo  diversos  procedimientos  para  ello.  Se  per¬ 
sigue  encontrar  para  la  posición  final  un  ND  que  ex¬ 
prese  el  valor  radiométrico  más  fiel  con  el  original. 
Los  procedimientos  más  usuales  son  el  del  vecino 
más  próximo ,  la  interpolación  bilineal,  y  la  convo- 
lución  cúbica. 


En  los  párrafos  siguientes  se  utilizará  la  denomina¬ 
ción  de  imagen  transformada  para  designar  a  la  re¬ 
sultante  de  la  transformación  geométrica  mediante 
las  funciones  de  ajuste,  e  imagen  corregida  a  la  re¬ 
sultante  del  remuestreo  de  la  transformada. 

Vecino  más  próximo 

El  método  del  vecino  más  próximo  (nearest 
neighbour)  consiste  en  asignar  a  cada  celda  el  ND 
correspondiente  al  de  la  posición  entera  más  cercana. 
Algebraicamente  se  puede  expresar: 


Ajc<0,5 

Ay  <  0, 5 

f{x+l,y) 

Av>0,5 

Ay  <0,5 

f(x,y+ 1) 

Ajc<0,5 

Ay  >0,5 

f(x+l,y+l) 

A*>0,5 

Ay  >0,5 

f(m,n)  = 


siendo: 

m  =  x  +  Ax  0<A.v<l] 
n  =  y  +  Ay  0  <  Ay  <  1 J 


(29) 


(30) 


donde  f{x, y)  es  el  ND  de  la  imagen  transformada  en 
una  posición  intermedia  determinada  por  un  número 
real ),f(m,ri)  es  el  ND  de  la  celda  (m,n)  en  la  imagen 
corregida  y  Ax,A y  representa  la  fracción  de  celda  en¬ 
tre  las  posiciones  (x,y)  y  ( m,n ). 

Este  procedimiento  no  altera  los  ND  de  la  imagen 
original,  sino  que  solamente  traslada  ligeramente  la 
posición  de  la  celda.  Presenta  el  inconveniente  de  po¬ 
der  introducir  fracturas  en  el  trazado  de  los  rasgos 
lineales  de  la  imagen  ( efecto  escalera).  Sin  embargo, 
es  el  único  método  utilizable  en  la  corrección  de  imá- 


Remuestreo  por  el  vecino  más  próximo. 
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genes  con  información  cualitativa,  como  son  las  imá¬ 
genes  clasificadas,  debido  a  que  respeta  el  valor  de  la 
celda. 


Considerando  coeficientes  de  ponderación  a  para 
cada  una  de  las  cuatro  posiciones  más  cercanas  a 
( m,n ),  se  puede  escribir: 


Interpolación  bilineal 

La  interpolación  bilineal  asigna  a  la  celda  corregida 
la  media  ponderada  de  los  ND  correspondientes  a  las 
cuatro  posiciones  más  cercanas  de  la  imagen  trans¬ 
formada.  El  peso  asociado  a  cada  nivel  digital  es  pro¬ 
porcional  a  la  proximidad  entre  ellos  (1-Ax,  1-Ay), 
medida  entre  centros  de  celdas. 

En  lugar  de  calcular  distancias  cuadráticas,  que  re¬ 
quieren  mayor  número  de  cálculos  de  proceso  infor¬ 
mático,  se  realiza  una  interpolación  por  filas,  obte¬ 
niéndose  primeramente  dos  interpolantes  con  los  va¬ 
lores  provisionales  en  la  fila  anterior  y  posterior  a  la 
posición  final.  En  un  segundo  paso,  se  realiza  una 
nueva  interpolación  entre  ambos  interpolantes  para 
deducir  el  ND  que  se  atribuirá  a  la  posición  ( m,n ). 


/  (m,  n)  =  c,  /  (x,  y)  +  c2f  (x+\,  y)  + 

+c3  /  (x,  y+l)  +  c4f  (x+l,  y+l) 

siendo  los  factores  de  ponderación: 

cx  =  (l-Ax)  (1-Ay) 
c2=  Ax  (1-Ay) 
c3  =  (1-Ax)  Ay 
c4  =  Ax  Ay 


(33) 


(34) 


Este  método  produce  resultados  suavizados,  pues  las 
cuatro  posiciones  más  próximas  contribuyen  al  ND 
final.  Es  más  apropiado  cuando  se  trata  de  corregir 
imágenes  con  información  cuantitativa. 


Los  interpolantes  /i  e  h,  correspondientes  a  la  fila 
anterior  y  posterior  a  la  posición  buscada,  tiene  los 
siguientes  valores: 

/  =(l-Ax)/(x,y)  +  Ax/(x  +  l,y) 

I2  =  (l- Ax)/ (x,  y  +  l)  + Ax  /  (x  +  l,  y  +  1) 

Y  el  ND  final  de  la  celda  corregida  será: 

/  (m,n)  =  (1  -  Ay)  /  +  A yl2  = 

=  (i- Ax)(i- Ay)/  (x,  y)  +  Ax(l- Ay)  /  (x  +  l,  y)  + 
+  (l  —  Ax)  Ay  / (x,  y  + 1)  +  AxAy  / (x  + 1,  y  + 1) 

(32) 


Convolución  cúbica 

La  convolución  cúbica  es  un  interpolador  que  emplea 
polinomios  de  tercer  grado.  Se  involucran  las  dieci¬ 
séis  celdas  más  cercanas  al  punto  considerado.  Los 
ND  se  interpolan  linealmente  en  grupos  de  cuatro  lí¬ 
neas  de  cuatro  píxeles  cada  una  para  formar  cuatro 
interpolantes.  Posteriormente  se  realiza  otra  interpo¬ 
lación  lineal  entre  los  cuatro  valores  obtenidos  para 
asignar  el  resultante  a  la  celda  corregida. 

Los  cuatro  interpolantes  se  calculan  del  siguiente 
modo: 


/(.v,y)^A.v 

A  j/(x+l,y) 

Ay 

f(x,y+l) 

1 

/(x+l,y+l) 

Remuestreo  por  interpolación  bilineal. 
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f{m)=  (A*)3- (A*)2  f(x  +  2)~ 

(A.r)3-(Ax)"-(A*)  f(x  + 1)  + 
(Ax)3-2(Ajc)2+1  f  (x)  ~ 

(A.v)3-2(A.x:)~+(Ajc)  f{x  —  l) 

Y  el  ND  final 


f(m,n)=  (Ayf-(Ayf  f(m  +  2)- 


( Aj)3-  ( Ay  )2-  ( Aj)  f(m  + 1) + 


(Ay)3-2(Ay)2+l  f(m)  - 
(Ay)3-2(Ay)2+(Ay)l/(m-l) 


(36) 


Modelización  matemática 

Cuando  las  distorsiones  sufridas  por  la  imagen  son 
sistemáticas  es  posible  conocer  la  función  que  ha  ori¬ 
ginado  la  deformación  y,  en  consecuencia,  la  imagen 
puede  ser  corregida  aplicando  algunas  transformacio¬ 
nes  simples. 

Traslación 

u  =  x  +  a 
v  =  y  +  b 


(37) 
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Es  un  procedimiento  muy  respetuoso  con  los  valores 
del  entorno  de  la  posición  corregida,  si  bien  es  el  que 
mayor  capacidad  de  procesamiento  requiere. 

Si  bien,  para  su  mejor  comprensión,  los  procesos  ex¬ 
plicados  han  seguido  la  dirección  imagen  transfor¬ 
mada— ^-imagen  corregida,  metodológicamente  se 
opera  al  contrario:  al  escribir  una  posición  entera  de 
columna  y  línea  en  la  imagen  de  salida,  hay  que  bus¬ 
car  el  mejor  ND  o  la  forma  de  calcularlo  a  partir  de 
las  posiciones  más  próximas  a  la  posición  intermedia 
en  la  imagen  transformada 


Inclinación 

u  =  x  +  ay | 
v  =  y  J 


(38) 


Cambio  de  escala 

u  =ax 
v  =  by 


(39) 
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Perspectiva 

u  =  axy] 
v  =  y  f 


(40) 


Rotación 

u  =  eos  dx  +  sen  9  y 
v  =  -sen0^:  +  cos9y 


(41) 


Si  bien  las  distorsiones  introducidas  en  la  imagen 
como  consecuencia  de  la  oscilación  de  la  plataforma, 
al  ser  aleatorias  suelen  corregirse  recurriendo  a  las 
transformaciones  polinómicas,  algunas  otras  pueden 
corregirse  mediante  modelización  matemática  apli¬ 
cando  operadores  sencillos.  Por  ejemplo,  para  el 
efecto  panorámico  puede  aplicarse  la  transformación: 


u 

( h/ x)  arctg  ( x¡ h)  9 

X 

V 

9  1 

_y_ 

(42) 


donde  h  es  la  altura  del  satélite. 

La  inclinación  de  la  órbita  se  corrige  haciendo: 

u 
v 

siendo  9  el  ángulo  formado  por  la  traza  del  satélite  y 
el  paralelo  del  lugar. 

La  ventaja  de  trabajar  con  operadores  matriciales  es 
que  se  pueden  combinar  varias  correcciones  simple¬ 
mente  multiplicando  las  matrices  que  caracterizan  a 
cada  una.  Por  ejemplo,  la  eliminación  del  desplaza¬ 
miento  de  líneas  debido  a  la  rotación  terrestre  y  la 


eos  9  senQ 
-senQ  eos  9 


(41) 
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inclinación  de  la  órbita  se  corrigen: 


u 

1  -a 

eos  9 

senQ 

X 

V 

0  1 

-senQ 

COS0 

_y_ 

eos  9  +  asenQ 

senQ-  a  cos0 

X 

- senQ 

COS0 

_y_ 

(43) 
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Traslación  Inclinación 


□ 


Cambio  de  escala 


Perspectiva  Rotación 


Transformaciones  geométricas  simples. 


